clc
clear
namePool = [0:2:28 29 30:2:82 83 84:2:144 145 146:2:194];
f = [0.0202:0.0202:0.303 0.312 0.3232:0.0202:0.8485 0.868 0.8687:0.0202:1.4747 1.482 1.4949:0.0202:1.9798] * 1e12;



[Ip, Im, IT, IQ, IM] = FC_RadiativePower_singleU('.\SC\SC_',namePool, f);
YW_Plot_Fast_singleU({f / (1e12), Ip.Ip, f / (1e12), Im.Im, f / (1e12), IT.IT, f / (1e12), IQ.IQ, f / (1e12), IM.IM}, {'-o', '-d', '-v', '-s', '-^b'}, ["\fontsize{28}\fontname{Helvetica}\itP", '\fontsize{28}\fontname{Helvetica}\itM', '\fontsize{28}\fontname{Helvetica}\itT', '\fontsize{28}\fontname{Helvetica}\itQ_e', '\fontsize{28}\fontname{Helvetica}\itQ_m'], ["Frequency (THz)", 'Radiation   power(a.u.)']);
% set(gca, 'yscale', 'log');
YW_Plot_Fast_singleU({f / (1e12), IT.ITx, f / (1e12), IT.ITy, f / (1e12), IT.ITz, f / (1e12), IT.IT}, {'-o', '-d', '-v', '-s'}, ["\fontname{Helvetica}\itTx", '\fontname{Helvetica}\itTy', '\fontname{Helvetica}\itTz', '\fontname{Helvetica}\itT'], ...
    ["Frequency (THz)", 'Radiation   power (a.u.)']);
set(gca, 'yscale', 'log');
